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Abstract 

In a recent paper [7] we have treated the future non-linear stability for reflection symmetric 
solutions of the Einstein- Vlasov system of Bianchi types II and VIo- We have been able now to 
remove the reflection symmetry assumption, thus treating the non-diagonal case. Apart from 
the increasing complexity the methods have been essentially the same as in the diagonal case, 
■^^i showing that they are thus quite powerful. Here the challenge was to put the equations in a 

( — ^ . form that permits the use of the previous results. We are able to conclude that after a possible 

■ basis change the future of the non-diagonal spacetimes in consideration is asymptotically 

diagonal. 



1 The Einstein- Vlasov system 



A cosmological model represents a universe at a certain averaging scale. It is described via a 
Lorentzian metric g a p (we will use signature — h + +) on a manifold M and a family of fundamental 
observers. The metric is assumed to be time-orientable, which means that at each point of M the 
two halves of the light cone can be labelled past and future in a way which varies continuously 
from point to point. This enables to distinguish between future-pointing and past-pointing timelike 
vectors. This is a physically reasonable assumption from both a macroscopic point of view e.g. 
\^ ' the increase of entropy and also from a microscopic point of view e.g. the kaon decay. One has 

also to specify the matter model and this we will do in the following. The interaction between the 
geometry and the matter is described by the Einstein field equations (we use geometrized units, 
£S| . i.e. the gravitational constant G and the speed of light in vacuum c are set equal to one): 

G a fi = 8-!rT a j3 

where G a p is the Einstein tensor and T Q( g is the energy-momentum tensor. For the matter model 
we will take the point of view of kinetic theory [12] . The sign conventions of [10] and the Einstein 
summation convention that repeated indices are to be summed over are used. Latin indices run 
from one to three and Greek ones from zero to three. 

We will consider from now on that all the particles have equal mass to. We will choose units 
such that to = 1 which means that a distinction between velocities and momenta is not necessary. 
The collection of particles (galaxies or clusters of galaxies) will be described (statistically) by a 
non-negative real valued distribution function f(x a ,p a ) on the mass shell. This function represents 
the density of particles at a given spacetime point with given four-momentum. Using the geodesic 
equations the restriction of the Liouville operator to the mass shell has the following form 

r) 8 

where are the components of the metric connection. We will consider the collisionless case 
which is described via the Vlasov equation: 

L(f)=0 
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The unknowns of our system are a 4-manifold M, a Lorentz metric g a p on this manifold and 
the distribution function / on the mass shell defined by the metric. We have the Vlasov equation 
defined by the metric for the distribution function and the Einstein equations. It remains to define 
the energy-momentum tensor T a p in terms of the distribution and the metric. Before that we need 
a Lorentz invariant volume element on the mass shell. A point of a the tangent space has the 
volume element \g^\^ dp°dp 1 dp 2 dp 3 (g^ is the determinant of the spacetime metric) which is 
Lorentz invariant. Now considering p° as a dependent variable the induced (Riemannian) volume 
of the mass shell considered as a hypersurface in the tangent space at that point is 

w = 2H(p a )6(p aP a + m 2 )\g^\Up°dp 1 dp 2 dp 3 

where S is the Dirac distribution function and H(p a ) is defined to be one if p a is future directed 
and zero otherwise. We can write this explicitly as: 

w = \p \- 1 \g^\idp 1 dp 2 dp 3 

Now we define the energy momentum tensor as follows: 

T a (3 = J f(x a ,p a )p a PfSTD 

One can show that T a p is divergence-free and thus it is compatible with the Einstein equations. 
For collisionless matter all the energy conditions hold. The Vlasov equation in a fixed spacetime 
can be solved by the method of characteristics: 

dXa _ pa. dpa _ r a p p pl 

Let X a (s,x a ,p a ), P a (s,x a ,p a ) be the unique solution of that equation with initial conditions 
X a (t, x a ,p a ) — x a and P a (t,x a ,p a ) = p a . Then the solution of the Vlasov equation can be 
written as: 



f(x a ,p a ) = f o (X a (0,x a ,p a ),P a (0,x a ,p a )) 

where /o is the restriction of / to the hypersurface t = 0. It follows that if /o is bounded the 
same is true for /. We will assume that / has compact support in momentum space for each 
fixed t. This property holds if the initial datum /q has compact support and if each hypersurface 
t = to is a Cauchy hypersurface. Before coming to our symmetry assumption we want to briefly 
introduce the initial value problem for the Einstein- Vlasov system. In general the initial data 
for the Einstein-matter equations consist of a metric g a b on the initial hypersurface, the second 
fundamental form k a b on that hypersurface and some matter data. Thus we have a Riemannian 
metric g a b, a symmetric tensor k a b and some matter fields defined on an abstract 3-dimensional 
manifold S. 

Solving the initial value problem means embedding S into a 4- dimensional M on which are 
defined a Lorentzian metric g a /3 and matter fields such that g a b and k a b ar e the pullbacks to S of 
the induced metric and second fundamental form of the image of the embedding of S while / is 
the pullback of the matter fields. Finally g a @ and / have to satisfy the Einstein-matter equations. 

For the Einstein- Vlasov system it has been shown that given an initial data set there exists 
a corresponding solution of the Einstein- Vlasov system and that this solution is locally unique 
up to diffeomorphism. The extension to a global theorem has not been achieved yet. However 
if one assumes that the initial data have certain symmetry, this symmetry is inherited by the 
corresponding solutions. In particular for the case we will deal with, i.e. expanding Bianchi 
models (except type IX) coupled to dust or to collisionless matter the spacetime is future complete 
(theorem 2.1 of 0). 
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2 Bianchi spacetimes 



The basis for the classification of homogeneous spacetimes is the work of Bianchi which was 
introduced to cosmology by Taub. Here we will use the modern terminology and we define Bianchi 
spacetimes as follows: 

Definition 1. A Bianchi spacetime is defined to be a spatially homogeneous spacetime whose isom- 
etry group possesses a three-dimensional subgroup G that acts simply transitively on the spacelike 
orbits. 

Our results concern only a special class of the Bianchi spacetimes, namely that of class A. 

Definition 2. A Bianchi A spacetime is a Bianchi spacetime whose three-dimensional Lie algebra 
has traceless structure constants, i.e. C% = 0. 

We will study II and VL> For Bianchi II the only non- vanishing structure constants are: 

Q23 = 1 = — C32 (1) 

and in the case of Bianchi VIo these are: 

C 31 = 1 = — C 13 , C 21 = 1 = — C 12 (2) 

We will use the metric approach. If W a denote the 1-forms dual to the frame vectors E a the 
metric of a Bianchi spacetime takes the form: 

4 g= -dt 2 + g ab (t)W a W b (3) 

where g a b (and all other tensors) on G will be described in terms of the frame components of a 
left invariant frame. A dot above a letter will denote a derivative with respect to the cosmological 
time t. We will use the 3+1 decomposition of the Einstein equations as made in [10]. Comparing 
our metric (J3J) with (2.28) of [TO] we have that a — 1 and (3 a — which means that the lapse 
function is the identity and the shift vector vanishes. There the abstract index notation is used. 
We can interpret the quantities as being frame components. There are different projections of the 
energy momentum tensor which are important 

p = T oa 
7 - T° 

J a — -L a 
S a b = Tab 

where p is the energy density and j a is the matter current. 
The second fundamental form k a b can be expressed as: 

g a b = -2fc a b. (4) 

The Einstein equations: 

kab = Rob + k kab - 2k ac k - 8n(S a b - ^9abS) - 4:npg a b (5) 

where we have used the notations S = g ab S a b, k = g ab k a b, and R a b is the Ricci tensor of the three- 
dimensional metric. The evolution equation for the mixed version of the second fundamental form 
is (2.35) of [TO]: 

k% = R a b +kkl- 87rS h a + 4tt<5? (S-p) (6) 

From the constraint equations since k only depends on the time variable we have that: 

R-ka b k ab + k 2 = 16ttp (7) 
V a fc a6 = 8iTj b (8) 
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where R is the Ricci scalar curvature. 

Another useful relation concerns the determinant g of the induced metric ((2.30) of |10]): 

|(log.g) = -2fc (9) 

Taking the trace of ©: 

k = R + fc 2 + 4tt5 - 12np (10) 
With ([7|) one can eliminate the energy density and (|T0|) reads: 

k=^(k 2 + R + 3k ab k ab ) + ATrS (11) 

Finally if one substitutes for the Ricci scalar with (Q: 

k = k ab k ab + 4^(5 + p) (12) 

Now with the 3+1 formulation our initial data are (gij(to), kij(to), f(to)), i.e. a Riemannian metric, 
a second fundamental form and the distribution function of the Vlasov equation, respectively, on 
a three-dimensional manifold S(to). This is the initial data set at t = to for the Einstein- Vlasov 
system. 

We assume that k < for all time following [8] (see comments below lemma 2.2 of [8j). This 
enables us to set without loss of generality to = —2/k(to). The reason for this choice will become 
clear later and is of technical nature. 

We will now introduce several new variables in order to use the ones which are common in 
Bianchi cosmologies and to be able to compare results. We can decompose the second fundamental 
form introducing a ab as the trace-free part: 

k a b = o ah - Hg ab (13) 



k ab k ab = <r ab a ab + 3H 2 (14) 

Using the Hubble parameter: 



we define: 



and 



Thus 





H = 




yb _ 








£_ = 


2E+ 


yl 


s? 




si 


si 



~(S5 + Ei) (16) 
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2 ^(SI-Si) (17) 



^3 
^3 



-E+ + \/3E_ 
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The reason for using the variables S + and £_ is that the diagonal case has been very important 
to understand the non-diagonal case. Define also: 

Q = 8np/3H 2 (18) 

H , \ 

d i- H «> 
The time variable r is dimensionless and sometimes very useful. From ([7]) we obtain the constraint 
equation: 

1 



ff3 -(iJ-<™0 = fi-i 

and from (|TT|) the evolution equation for the Hubble variable: 

*0O-i + s<£ + + (2D 

Combining the last two equations with ^ we obtain the evolution equations for £_ and S + : 

£+ = g[ 2fl ~ 3 6 ( g + fl|) - S + (3 + ^) + ^(3Sg + 35 3 3 - 25)] (22) 

Since we use a left-invariant frame / will not depend on x a and the Vlasov equation takes the 
form: 

p dt p dp*- 

It turns out that the equation simplifies if we express / in terms of Pi instead of p % what we can 
do due to the mass shell relation: 

P -_rWV— = o 

Because of our special choice of frame the metric has the simple form . Due to the fact that we 
are contracting and the antisymmetry of the structure constants we finally arrive at: 

g + ltf-'Cllft-g-O (24) 
From it is also possible to define the characteristic curve V a : 

^ = (V a )- l C d ba V b V d (25) 
for each Vi(t) = Vi given i. Note that if we define: 

V = gVViVj (26) 
due to the antisymmetry of the structure constants we have with (|25D : 

£ - i^ v ' v ' (27 > 

Let us also write down the components of the energy momentum tensor in our frame: 

Too = / f(t,p a )p°V9dp 1 d P 2 d P 3 (28) 

ley = -J f(t, p^PjVfjdp'dp^p 3 (29) 

Tij = I f{t 1 p a ) Pt p J (p )- 1 ^dp 1 dp 2 dp :i (30) 



5 



3 The asymptotics of Bianchi II and VIq 



Before coming to the non-diagonal case we have a look at the tilted fluid models, since they are 
non-diagonal as well and they may help us to understand the non-diagonal case with collisionless 
matter. For the tilted Bianchi II we use the corresponding equations of [5] and for Bianchi VIo 
the equations of [3], in both cases with 7 = 1. We will not go into the details for this we refer 
to the mentioned work. The point is that looking at the linearization we see that the variables 
which did not appear in the diagonal case have decay rates which are between the ones considered 
previously. This is a good sign. Also in [2] the stability of the Ellis-MacCallum solution, in fact 
the stability of the Collins solution, was already considered within the Einstein-Euler system. 

3.1 Equations of the non-diagonal case 

Using (|15p we arrive with ^ for a 7^ b to: 

^ = *[|t-^(3 + Ji)-^]; a*b 
which together with ([22 ]) -P5 ]> . i.e. 



describe the evolution of The expression for the Ricci tensor is: 

% = -\clAC% + 9 lm g kn C^ - \c^C p ql g jm g ipg ^g ln (31) 

and 

4 = R lb g b] = -\c l ki g hj [Ci + gi m g kn CZ) - \c'LC p ql g iP 9 kq 9 ln (32) 

We will now derive some expression concerning the derivative of (|31| : 

Rij = CLCnj'(^im.9 fcTl — gimk kn ) + 

2^rikC'qi(kj m gipg kq g ln + gj m ki P g kq g ln — gjm9ipk kq g ln — gj m gipg kq k ln ) 

Thus: 

g^Rij - gi r C l kl C™(k lm g kn -g lm k kn ) + 

\c P ql [C™ k k r m g lp g kq g ln + C r nk (k w g kq g ln - g lp (k kq g ln + g kq k ln ))} 

For r = i and relabelling the m with i for the terms with the prefactor ^: 

.</'"/>-, - 9 3l C l kl Cl) (k lm g kn - g lm k kn ) + \c p ql & nk [2k ip g kq g ln ~ g w (k kq g ln + g kq k ln ))) 
Rearranging terms: 

.</"/>,, = CLC™ (k lm g kn g^ - g lm k kn 9 n + C p ql C* nk {h p g kq g ln - g lp k kq g Ln }. 
We see that the first with the third and the second with the fourth term cancel each other, hence: 

gTRu = (33) 
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The evolution equation for the Ricci scalar due to ([55)) is: 

R = 2R)ki = 2H(-R + Rfii) 

Define 



R 3 

N 1 = -2- 

• H 2 



The derivative of this expression is: 



N{ = 9 -j^ + 2H(NfS J p - (1 + —)Nl) 

Consider the quantity N = R/H 2 . Its evolution equation is: 

N = 2H[qN + Np{\ (34) 

3.2 Curvature expressions 

For bookkeeping reasons we define the following quantities where we use from now on g for the 
determinant of the metric. 



.4 


= g 22 g 33 


- (3 23 ) 2 = 


5ii. 


B = 


5 13 5 23 - 


- 9 12 9 33 
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9 








9 


C 




- 9 U 9 22 = 


513. 


D = 


5 12 5 13 


-5 U 5 23 


_ 523 
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E 


= <? n 5 33 


(9 13 ) 2 = 


522 . 
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F = 


5 U 5 22 


- (5 12 ) 2 = 


_ 533 
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Let us denote the quantities divided by H 2 with small letters, i.e. a = -j^. 

3.2.1 Curvature expressions for Bianchi II 

Using (JjJJJ) for Bianchi II: 



Rl = \9iACUg 23 9 23 ~ 9 22 g 31 ) + C4(5 23 5 3J - 9 33 9 23 )] + ^9..(%A 



We obtain: 



and as in the diagonal case: 



r> 1 4 l(5ll) 2 

R = -2 mA = -2 — 



pi d r?2 r?3 

it^ — — it — — if.2 — — ii3 

c>2 r?3 d3 r>2 ri 

it^ — it^ — it 2 — it.3 — U 



However in the non-diagonal case we have: 



, 512 

'511 



Thus 



3 5ii 



iV = -2H[(l + -fj + ^+)N - Wu] 



7 



where Wn = N 2 T,f + iVgEf . In order to calculate the derivative of N 2 we need the following 
expression: 

i?^(-2— ) = 2H[2Y.\R + (3E + + v^-)^ - ^-{(R l 2 f^\ + ii^E?)] 
at gn zii 



Hence: 



Nl = H[ANV\ - 2(S+ + 1 - \/3£_ + -^)7V 2 1 + W%] 
Nl = H[4NVl - 2(£+ + 1 + V3E_ + -^)7V 3 1 + W^} 



where 



Wl = -2S^ 2 1 + AriJV-^S^ 1 + Y>lNl) 

3.2.2 Curvature expressions for Bianchi VIo 
With (|S2D we obtain: 

-2^ = ^(Cf, + Cl) + g a (C{ z E - C[ 2 D) + g l3 {-C{ 3 D + C{ 2 F) 
•//:;:< f; ! ,; ( r',y :: + g l > g 12 ) + (%■ .fUf' + g^g 22 )} 

■<n,Ciuyur g 2j g n ) + C 2 3l (g^g 23 - 5 2 V 3 ) + (%uy>^ - g^g 11 ) + CUg^g 23 g 3j g 21 )} 

In particular: 

R = -\[{V^2~E + ^fg^F) 2 - Ag 23 D] = -^[(<?22 + gss) 2 - 4 5 | 3 ] 

Rl = U922E - g 33 F) = ^[(g 22 ) 2 - (g 33 ) 2 ] 
2 2g 



and like in the diagonal case: 



However we have 



R — R\ 

r>2 d3 

K 2 — -H 3 
R\ =Rl = 



Nl = -Nl = g 23 (f - e) = -N 23 (N 3 + N 3 ) 

^1 = ~ 2 Jp + 9i2( e - /) = ^12(^2 - N 3 ) - 2N 13 N 23 

N i = - 2 Jj2 + 9i3(f - e) = N 13 (N 2 - N 3 ) - 2N 12 N 23 



where N t j is defined as 



Nij *L 



V9H 
and 

N 2 = N 22 
N 3 = ~N 33 
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which means that N 2 = R 2 /H 2 : 

Ni = lm) 2 -(N 3 f) 

Recalling that 




we can compute the derivatives of Nij using the following formula 

Nij = H[qN l3 - 2E?JVy] 

Hence 





N 12 = 


H[(q- 


4S+)iVi2 


- 2Y? X N 2 ~ 


2ZlN 23 ] 


(35) 




Nl3 = 


H[(q- 


4S+)iVi 3 




h 2S?A^ 3 ] 


(36) 


^23 


= tf[(2E + -\ 


-2V3Y,. 


- + q)N 23 


+ 2E|iV 3 - 




(37) 


N 2 


= tf[(2E + 4 


-2V3T,. 


- + q)N 2 - 


- TZ\N 12 - 


2E^ 23 ] 


(38) 


N 3 


= # [(2E+ - 


- 


- + q)N 3 - 


V 2ElN 13 + 


2E2jV 23 ] 


(39) 



From ([3~i)l we obtain 

AT = 2if[(2E+ + q)N - 2V3E_iVf + iVf E| + N^j + N^l + A^E^] (40) 
The evolution equation for iV|: 

N% = F[2(2E+ + g)JV| + 2V3E_((7V 3 ) 2 + (A^) 2 ) - 2(E^ 12 iV 2 + ^ 3 N 13 N 3 + Y* 3 2 N 23 N 2 + T,jN 23 N 3 )) (41) 

3.3 The non-diagonal asymptotics of Bianchi II and VI 

We will now discuss the asymptotics of the non-diagonal case. The structure of the analysis is 
very similar to the diagonal case. We start with a bootstrap argument and end with applying 
Arzela-Ascoli. Next we will collect the bootstrap assumptions. The pref actors denoted by A and 
some index are small constants. 

3.3.1 Bootstrap assumptions for Bianchi II 
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3.3.2 Bootstrap assumptions for Bianchi VIq 



1 


< 




-r 


B 




< 


A_(H 


-*r 


3 
S 


9 

l^ + ^i 

1 8' 


< 


4:1 (1- 


1- 1)~ 


3 
8 


l ?vr2 1 
1^2 1 


< 






3 
S 


1 AT 1 


< 


<~ 1 






1 a j 1 


< 


C2 






|Aaa| 


< 


4 23 






P 


< 


A m (l- 




1 

3 


|e§I 


< 


A 23 (l- 




3 
4 


l^| 


< 


A 32 (l- 




3 
4 


l^| 


< 


Ai 2 (l- 




3 
4 


l^| 


< 


Ai 3 (l- 




3 
I 




< 










< 


c 4 







3.3.3 Mean curvature 

Concerning the estimate of H there is no difference with respect to the diagonal case. The reason 
is that the estimate of D 

„ 1 / 3 „ h . AttS 
D = T2^ N+ IP^ + 3IP 

is the same. Thus as in the diagonal case it follows from fllTTl) that 

d t (H- 1 ) = ^ + 0(et-i) 
and following the steps made for the diagonal case we arrive at: 

H= 2 -t- 1 {l + 0(et~i)) 

will hold. 

3.3.4 Estimate of the metric and P 

For a matrix A its norm can be defined as: 

||4| = S vv{\Ax\/\x\ : x 56 0} 

Let B and C be n x n symmetric matrices with C positive definite. It is possible to define a 
relative norm by: 

||fi||o = sup{|Bx|/|Cx| : x 0} 

Clearly: 
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It also true that: 



\\B\\ < \\B\\ C \\C\\ 



\\B\\ C < V^(.C~ 1 BC- 1 B) (42) 

This can be shown as follows. Consider the common eigenbasis hi of B and C . Then there exist 
Ui such that Bbi = uiCbi for each i. Then (|42j) is equivalent to the statement that the maximum 
modulus of any is smaller than Y^a 2 . Using (|42p we obtain in the sense of quadratic forms: 

<r ab < (<i cd a cd )ig ab (43) 

Define 

g ab = th ab 

Then 

-|(i- V b ) = ^ 7_1 9 ab (-7 + J) + 2 ^ 7+ * K b - Hg ab ) 

where we have introduced for technical reasons a small positive parameter 7. Using now the 
inequality (|4*5)) 

J t (t~ 7 9 ab ) < t^- l g ab [-l + l q + 2tH((H- 2 a cd a cd )^ - 1)] (44) 
Using the equation f4"4")) and the estimate of H 

f t (t-^g ab ) < r-^i-T + v - + |(i + o^))^- 2 ^"^ - 1)] 

We obtain decay for the metric (in the sense of quadratic forms) provided that (H~ 2 <r c d<7 cd )^ < 1. 
This holds for Bianchi II and VIo with for instance | = 0.4. Thus we have 

9 ab < t-hlg ab (t ) 

This implies that the components of the metric are also bounded by some constant C(to) which 
depends on the terms of g ab (to). Consider now 
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2H{Y b a - S b a )g a f 



Since the metric components are bounded the non-diagonal terms will contribute only with an e. 
Thus we have for every component g lJ (no summation over the indices in the following equation): 

ff« = 2H(^ i - 1 + e)gV < 2ff (max(£*) - 1 + e ) 5 « = 2H{~ + e)g^ 
Using now the estimate of H 

g ij <t- 1 (-l + e)g ij (45) 

One can conclude that 

Wg-'W <or 1+£ ) 
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From (gSj 

V = g bf V b V f < r x (-i + e)V 

which means that 

V = 0(^ 1+£ ) 

which gives us the same decay for P as in the diagonal case: 

P = 0(t"3 +e ) 

3.3.5 Closing the bootstrap argument for Bianchi II 

It follows immediately by the same arguments as in the diagonal case: 

£_ = 0(t- 1+e ) 
E 2 = 0(t- 1+e ) 

sf = o(r 1+e ) 
s| = o(r 1+£ ) 



E3 = or 1+e ) 



Defining (iVi) 2 = — 2iV we arrive at: 

£+ = H[@f- £+(3 + A) + _g_( 3 <?f + 3S| - 25)] 
E_ = g [-(3 + 4)E_ + ^^l] 

Since 2-^ L decays like t~ 1+e we see that we can apply the same arguments as in the diagonal case 
to obtain an improvement of the bootstrap assumptions: 

The system which remains using the time variable r is the following: 

(Ei)' = E i( ff -2)+JV£-5g» 

(E5)' = Ei(g-2) + JVj-^ 

(JVa 1 )' = -2(7V 1 ) 2 Ei - 2(E+ - 9 - v^E-)^ 1 + W\ 
= -2(JVi) 2 Ej - 2(£+ - q + V3X_)N% + W£ 

Let us focus on the Y,\ — iV^-system. Using the estimates obtained we arrive at: 
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Let us go to the basis of eigenvectors of the linear system via the linear transformation 
Thus we arrive at 

§)'-(-<? sxd+^-^C 

Using the bootstrap assumptions for £ 2 and iV 2 we have an assumption for £ 2 . By the usual 
contradiction argument we arrive at 

Integrating the equation for iV 2 we arrive at 

Nl = Nl{r ) + O{e) 
Going back to the variables £ 2 and N% via 

£|\ _ 1 (A 2\ (t\ 



Nl) 3 V3 3 J \Nl 

Ej(r) = [2Si(r )-liV 2 1 (ro)]e(-i + ^ + ^ 2 1 (r )-Si(To) + 0(6) 
Ni(r) = [^i(r )-iV 2 1 (r )]e(-i + ^ + 27V 2 1 (ro)-^i(r ) + 0( e ) 
Changing back to the time variable t: 

Y.\{t) = C(to)[2I]i(to)-^iV 2 1 (^)]t-^ +£ + ^iV 2 1 (io)-Si(t ) + 0( e ) 
N}(t) = C(to)[|sJ(t ) - ^(to)]*-^ 8 + 2tf 2 1 (t ) - |eS(*o) + 0(e) 

where C is a constant, in particular C(to) = t£ e * T ° . The only term which could prevent us from 
improving the estimates is the e coming from the bootstrap assumptions of £ 2 , but note that it 
comes in combination with £^ as a product of both, thus the last term O(e) on the right hand 
side of the last two equations does not prevent us from improving our estimates. Thus if we wait 
long time enough and choose iV 2 (io) and £ 2 (io) small enough we will have an improvement for 
iV 2 and £ 2 since we can choose them independently and smaller then Ayi and A c ±2- There is no 
difference in the procedure for N% and £3. 

3.3.6 Arzela-Ascoli for Bianchi II 

Since all estimates have been improved we can apply Arzela-Ascoli and we arrive for £ 2 and iV 2 
to: 

£*(t = oo) = Itf^to) - £j(tb) 
N%(t = oc) = 2N 2 1 (t Q )-^ 1 2 (t ) 

Consider now the following transformation of the basis vector 

£1 = ei 

e 2 = e 2 + aei 

§3 = e 3 + be\ 
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It preserves the Lie-algebra, i.e. the Bianchi type. The following relation holds between the 
variables E5, an d E3 in the different basis: 

'El Ef E?\ /l 0\ /E} \ / 1 0\ / El 

E 2 E| Ef = a 1 J I E^ -a 1 = E 2 + a(Ej - E|) E| 
,Ei S| E|/ V & V V S 3 Ei/ V-6 l/ V S 3 + fe(S{-E|) E|, 




^3/ 

We see that choosing a = — ^2(00) and b = — |Eg(oo) the transformed variables Ej, E3 are 
zero asymptotically. By direct calculation one can see that the same is true for the transformed 
variables TVj and N^. Thus we obtain the same asymptotics as in the diagonal case and we can 
conclude: 

Theorem 1. Consider any C°° solution of the Einstein- Vlasov system with Bianchi II symmetry 
andwithC 00 initial data. Assume that |E + (i ) - 1|, |E_(< )|, |E^(t )|, [53^(t )|, |E§(t )|, |El(i )|, 




3 
4 

times, after possibly a basis change, the following estimates hold. 



T,l(t )\, |E?(f )|, |JVi(* ) ~ ||, \N%(t )\, \N£(t j\ and P{t ) are sufficiently small. Then at late 



H(t) = l t -\l + 0(t-i)) 

E + -i = 0(H) 

E_ = Oit- 1 ) 

E 2 = 0(H) 

E 3 = 0(H) 

E* = 0^) 

Ei = Oit' 1 ) 
Ef = Oit' 1 ) 

Ni~l = Oir--) 

N 1 = 0(H) 
Nl = 0(H) 
Pit) = 0(H) 

3.3.7 Closing the bootstrap argument of Bianchi VIo 

It follows immediately by the same arguments as in the diagonal case: 

) (46) 
) (47) 
) (48) 

Now consider the E2-/V23 system. Using the fact that iV| = —N 23 iN 3 + N 2 ) we obtain 

Si = H[-N 23 iN 3 + N 2 ) St(3 + Jy) - ^gt] 

AT 23 = H[i2E+ + 2y/SS- + q)N 23 + 2E%N 3 - 2E 2 iV 13 ] 



IE2I 


= 0(t" 


-l+e 


IS3I 


= o(r 


-1+e 


£2+£§| 


= 0(t" 


-l+e 
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Using the bootstrap assumptions, the estimates obtained and the variable r: 



y3 

^2 



23 



£1 £2 
£3 £i 



V3 
-2 



N; 



2:; 



0(ee 



where ei, e2 and £3 have the following origin. The quantity t\ is determined essentially by the 
error in N and E + and note that E| comes in combination with E§, thus this term can be chosen 
as small as we want. The quantity £2 comes from N2 + N3 and can be determined by the error 
of JV| and finally the quantity £3 which comes from N3 depends on the error of N, N$ and 7Vf 3 . 



Note in the last term that the quantity is squared, thus it is negligible, 
linearization and going to the eigenbasis via 



Having a look at the 



-2 



N- 



23 



1 

-1 1 



V3 
-2 



23 



we come to the system 



y3 

—2 

.N 23 

From which follows 



+ £1 + £2 
£3 - £2 



£2 
£1 - £2 



N- 



23 



0(ee 



N 23 = TV 23 (t ) + 0( e ) 



and going back 

E|: = S|(ro)e(-l+^ 

N 23 = iV 23 (T )-E3(T ) + O(£) 

We see that we have improved N23, ^2 an d with that also E§ 

Using these estimates and the bootstrap assumptions let us focus now on the following system: 



E_ 



= H 



HI 



s+(9-2) + o(r 1+f )] 

AT 2 

-^| + ( (Z -2)E_+Or 1+e )] 



N = ff[2(2E + + q)N - 4\/3E_iV 2 2 + +0(t^ 1+€ )} 



N. 



i>2 



= H [2(2E+ + q)Nl + (-V3 + 0(e))E_ + 0(<- 1+e )] 



where in the last equation (iV 2 ) 2 + (-/V3) 2 was estimated with iV|, TV and 7V 2 3. The 0(£)-term will 
not play a role since it can be absorbed in the e of the estimate. Let us look at the linearization 
using the variables E + = E + + ±, E_ = E_, N = N ■ 











N 




E_ 




\N$) 


V 



21 

II 

16 



►4 

16 









_3 



£ TV 2 

8' JV 2 






TVj and and the variable r 



\ 






N 




E_ 


/ 


U 2 V 



The eigenvalues are 



<Vi/ 
A 



3/4 



-7 ±7^ 
4 4 

3 3. 
— ± —1 

4 4 
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These eigenvalues are the same which appeared in the reflection symmetric case. Using the same 
arguments we arrive at 



\Z + + \\<A + (l + ty^ 
|S_| < A-(l + t)-i +e 
\N+l\ <A cl (l + t)-^+ e 



\Ni\ <A c2 (l + ty 



Finally, only Y^(tq), Ef (tq), Ni2(tq) and Niz(tq) have to be improved. Let us look at the Ef , N12 
system. There is no difference between this system and the Ef-^13 system. 

0=(:t- f:iO(£W>C 



Ax 
A 2 



= -(2ei - V2V2e 2 e 3 - 3e 2 + 3e 3 ) 
= i(2 ei + \/2V2e 2 e3-3e 2 + 3e3) 



Now choosing the error of JV bigger than E ab E a {, t\ will be negative. e 2 and €3 can be chosen in 
such a way that the square root of the term is positive but in total smaller than e\ , such that we 
have two small and different eigenvalues. 



3.3.8 Arzela-Ascoli for Bianchi VI 

For Bianchi VLj we can apply Arzela-Ascoli as well. We see that N23 will be zero and E^(to), 
Ef(ro), A r i 2 (ro) and Ni^{tq) will be come constants. This time we can make the following basis 
change which preserves the Lie-algebra to obtain that the mentioned variables tend to zero: 

ei + ae 2 + be 3 
('2 
('?, 



ei = 
e 2 = 
63 = 



We can conclude 



Theorem 2. Consider any C°° solution of the Einstein- Vlasov system with Bianchi VIq symmetry 
and with C°° initial data. Assume that |E + (< ) + j\, |E-(*o)|, E|(t ), |E£(t )|, |E|(t )|, |Ef (i )|, 
|£?(i )|, |S?(t )|, \N(t ) + ||, \Ni(t )\, |E?(t )|, |JVia(to)|,|JVi 3 (to)| and P(t ) are sufficiently 
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small. Then at late times, after possibly a basis change, the following estimates hold: 

0(t-i)) 



tl[t) 


YVl 


1 

4 


= o(t-i) 


£_ 


= 0(t~i) 




= Oir 1 ) 


=5 


= oit- 1 ) 


^3 


= o^- 1 ) 


V 3 

-2 


— ^ \~ ) 


^1 


— n(+—i\ 


=? 


= 0(t~i) 


iVl2 


= 0(t~i) 


JVl3 


= 0(t~i) 


V 23 


= 0(t~i) 


TVf 


= o(t~i) 


9 

N+ l 


= 0(t~i) 


P(t) 


= 0(*~*) 



4 Conclusions 

As mentioned in the abstract the challenge here was to put the equations in a form such that the 
results of the diagonal case can be used. This can be seen especially in the curvature variables. 

R: 



. For Bianchi VIo we had to 
The notation might be a little bit 



For Bianchi II it was sufficient to use the new variables Nj 
introduce in addition to that the new variables Na — 9 %- . 
confusing, but in both cases these variables have a connection to the curvature variables N±, N2 
and 7V3 of the diagonal case and this is the reason for the notation. In contrast to the diagonal 
case where the treatment of Bianchi II and VIo was almost identical, here the latter case was more 
difficult. One reason could be the obvious increase in complexity. In the Bianchi II case it was 
sufficient to deal with TV instead of TVi and look at the differences. In the case of Bianchi VIo 
TV had to be used to start the bootstrap argument. Then also TV| and TV23. This last variable 
made the correspondence to the diagonal case more difficult. As can been seen in the chapter 
where the bootstrap argument was closed for Bianchi VIo, we had to look more carefully on the 
dependence of the different e. Note also that we did not use exactly the linearization in our last 
improvement of the estimates. We would have obtained that zero is an multiple eigenvalue and we 
would have not obtained decay, but logarithmic growth. This would have been sufficient to close 
the bootstrap argument with corresponding suitable bootstrap assumptions, but there would exist 
difficulties to apply the Arzela-Ascoli theorem and to obtain that the non-diagonal components 
become constant. Another difference to the diagonal case is the use of a basis change in the end. 
In general the non-diagonal components will become constants and thus not relevant. However to 
obtain "diagonal" asymptotics a basis change will in general be necessary. 

It would be interesting to investigate whether the work on homogeneous Ricci solitons [3] 
can help to understand the similarities and differences between Bianchi II and VIo (in Thurstons 
classification Nil and Sol). 

We have discussed the future asymptotics of some Bianchi models, what about the higher types? 
The case of Bianchi VIIo will probably be quite different. For instance in [T3J it was discovered 
that the Bianchi VIIo spacetimes with a non-tilted fluid are not asymptotically self-similar in the 
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future and that some oscillations take place. It is shown that dynamics are dominated by the 
Weyl curvature. However for dust a bifurcation of the Weyl curvature takes place (theorem 2.4 
of |13| and comments below). For this reason it is likely to expect difficulties when applying our 
techniques to this case. Something similar, but even more complicated happens in the case of 
Bianchi VIII spacetimes with a non-tilted fluid [5]. 

What about inhomogeneous models? Some direction to generalize our results could be to 
analyze the Gowdy model which is the simplest inhomogeneous case. In |llj different links between 
Bianchi and (twisted) Gowdy spacetimes are considered, in particular for Bianchi I, II, VIo and 
VIIo- The analysis of perturbations is another interesting approach towards the understanding of 
inhomogeneous models (see PQ for recent developments). 
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